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SECTION  I 


INTRODUCTION 


A.  OBJECTIVES 

The  primary  objective  of  this  study  was  to  investigate  the 
gravitational  effects  in  an  elastic  medium.  More  specifically',  the 
objectives  of  this  study  were  to: (1)  derive  a  governing  equation  of  the 
motion  with  gravity  force  term,  (2)  develop  an  experimental  technique 
to  simulate  the  gravitational  effects,  (3)  determine  the  parameters  that 
govern  the  gravitational  effects  in  an  elastic  medium  using  a  modi  1 
technique  in  the  laboratory,  and  (4)  gain  experience  necessary  for 
detailed  model  studies  in  the  Phase  II  Program. 

B .  BACKGROUND 

Many  studies  have  been  performed  to  investigate  the  soil-structure 
interaction  problems  subjected  to  dynamic  loading  conditions. 
Research  work  simulating  this  problem  used  either  full-scale  or 
small-scale  laboratory  modeling  techniques.  ine  full-scale  technique 
with  the  necessary  instrumentation  (i.e.,  soil  stress  meters,  pore 
water  pressure  transducers,  settlement  gauges,  and  strain  gauges, 
etc.)  gave  the  best  results  for  estimating  prouoL^pc  behavior. 
However,  this  technique  has  shown  serious  major  drawbacks:  mainly, 
cost  and  time  of  construction  and  operation.  For  these  reasons, 
small-scale  modeling  has  been  used  as  an  alternate  testing  method. 

However,  the  use  of  small-scale  modeling  in  the  laboratory  has 
shown  severe  limitation  when  the  gravity  body  force  of  the  structure 


1 


itseli  is  the  principal  load  on  the  system.  This  is  because:  (1)  soil 
characteristics  are  nonlinear  and  overburden  dependent,  and  (2)  smaller 
stress  magnitudes.  The  stresses  in  a  small-scale  model  due  to  its  own 
weight  are  much  smaller  in  magnitude  than  those  in  the  corresponding 
prototype  system.  Because  of  this  limitation,  the  effects  of  the 
gravity  has  been  neglected  in  the  small-scale  model  technique  based  on 
the  assumption  that  gravitational  and  strain  rate  effects  are 
insignificant . 

In  recent  years,  however,  a  new  model  technique  (centrifuge  model 
technique)  has  been  introduced  and  applied  to  small-scale  models  to 
simulate  gravitational  effects.  Based  on  the  centrifuge  model  studies 
(References  4,  10  &  11),  it  has  been  demonstrated  that  the  gravity  scale 
is  an  important  parameter  in  small-scale  model  study.  The  purpose  of 
this  research  was  therefore,  to  determine  the  parameters  that  govern  the 
gravitational  effects  without  using  the  centrifuge  model  technique. 

C.  SCOPE 

This  contains  a  review  of  elastic  waves  in  a  solid  medium  without 
gravny  force;  a  derivation  of  a  governing  differential  equation  of  the 
motion  with  gravity  force;  a  description  of  developed  laboratory  testing 
facility;  and  preliminar>  model  rest  results  and  discussion.  Finally, 
conclusion  are  drawn  and  used  to  offer  recommendations  for  the  direction 
of  future  research. 


SECTION  II 


WAVE  PROPOGATION  IN  SOLIDS 

A.  INTRODUCTION 

Dynamic  load- generated  waves  can  be  divided  into  three  mail', 
categories:  compressive,  shear,  and  surface.  The  three  main  wace  types 

can  he  divided  into  two  varieties:  body  waves,  which  propogato  through 
the  bodv  of  soil,  and  surface  waves,  which  are  transmitted  along  a 
1  <«<.v .  The  most  important  surface  wave  is  Rayleigh.' s  wave.  Body 
v,.ve s  can  he  further  subdivided  into  compressive  or  longitudinal  waver;, 
and  distort  ior.al  or  shear  waves.  This  section  reviews  these  three  types 
o  i  w  a  v  t  •  . 


B.  EQUATION'  OF  MOTION  IN  AN  ELASTIC  MEDIUM  WITHOUT  GRAVITY  FORCE 

Figure  1  shows  an  isotropic  elastic  medium  whose  sides  measure  dx 
and  dv.  The  normal  stress  acting  on  the  planes  normal  to  the  x  and  y 
axes  are  a  and  a  ,  respectivel” .  The  shear  stresses  are  r  and  r 

x  >  icy  yx 

Let  u  and  v  be  the  displacements  in  the  direction  of  x  and  y, 
respectively.  Then  the  expressions  for  strains  and  rotations  of  the 
medium  are  : 


3u 

(  —  - 

X 

dx 


dv 

€  -  - 

y 

6y 


(1) 


(2) 


3 


cr 

y 


Figure  1.  Notations  for  Normal  and  Shear  stresses 
in  Cartesian  Coordinates  System 


Figure  2.  Derivation  of  the  Equations  of  Motion 
in  an  Elastic  Medium 
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where  , 


e  and  ( 


xv 


3u 

dv 

ru 

3y 


normal  strains  ii.  the  direction  of 
x  and  y,  respectively, 
shearing  strain,  and 
the  components  of  rotation  about  the 
x  and  v  axes 


(3) 


(4) 


The  generalized  Hooke's  law  is  used  for  the  relationship  between 
stresses  and  strains: 

1 

«  -  -  (o  -  v  o  )  .  (5) 

x  E  -  y 

1 

e  =  -  ( o  -  v  a  )  .  (6) 

y  y  x 

E 

T  G  ■  7  .  O) 

xy  xy 


where  E  is  Young's  Modulus,  v  is  Poisson's  Ratio,  and  G  is  shear 
modulus  defined  as  G  -  E  /  (2  (1  +  i/)). 

Equations  (5)  and  (6)  can  be  rewritten  in  terms  of  normal 
strains  as : 
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where , 


A-f  t-  2  G-e 

X  X 


a  —  A  ■  (  +  2G  •  e 

y  y 


(8) 

(9) 


A  =  »/-E/((l+i/)(l-i/))  .  (10) 

e  =  t  +  £  -  volumetric  strain. 

X  y 


To  obtain  the  equation  of  motion  for  an  elastic  medium,  stresses 
acting  on  an  element  with  sides  measuring  dx  and  dy  are  considered  as 
shown  in  Figure  2.  Four  separate  forces  will  be  acting  parallel  to  each 
axis,  and  if  one  considers  the  resultant  force  acting  in  the  x- 
direction : 


da 

X 

(-  a  +  (o  +  -  dx  ))dy 

X  X 

dx 


+  (-T  +  (r  + 

yx  yx 


dr 

yx 

-  dy ) )  dx 

dy 


which  simplifies  to 


(  — —  +  — ^  )  dx  dy 

dx  dy 


and  by  the  Newton's  second  law  of  motion,  neglecting  body  forces  such  as 

2  2 

gravity,  this  will  be  equal  to  (p  dx  dy  (3  u  /  dt  ))  where  p  is  density 
of  the  element,  so  that 


da 

X 


dx 


dr 

yx 


dy 


dZu 


P 


dt 


2 


(ID 


6 


S imilarly , 


do 

y 

+ 

dr 

xy 

a2v 

(12') 

ay 

<3x 

T 

Co 

rt 

to 

.  \  *•  / 

These 

equations 

of 

motion  will  hold, 

whatever  the  stress  -  strain 

behavior  of  the  medium.  In  order  to  solve  the  equations,  let  us 
substitute  Equations  (7)  and  (8)  into  the  Equation  (11) . 

It  yields 


or 


a‘  u  3  3  <3v  3u 

f  - — —  =  — —  (At  +  2G-e  )  +  G  -  (  -  +  -  ) 

X 

8t  3x  3y  ox  3y 


a2u 

3  € 

-  \  _l  n 

a2u 

(  - 

ax2 

a2v 

4- 

a2u 

a2u 

-f. 

P 

at2 

A.  T  VJ 

dx 

3x3y 

„  2 
ox 

ay2 

Since 


a2u 

_L 

a2v 

d  e 

.  2 

OX 

T 

dxdy 

dx 

the  equation  of  the  motion  becomes 

32u  3e 

p  -  -  (A  +  G)  -  +  G-VZu 

at2  ax 

where 


ax2  ay2 


Similarly , 


(13) 
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a2 
a  v 


3t 


9£ 

=  (A  +  G)  -  +  G  V  v 

dy 


(14) 


Now,  differentiating  Equations  (13)  and  (14)  with  respect  to  x 
and  y,  respectively,  and  adding 


„2 


a2'  a2' 

3  £  at 


3u  3v 


ad  u  3v 

p  ( -  (  -  +  -  ))  “  (A  +  G)(  -  +  -  )+  G  V  (  -  +  ) 


3t  dx  dy 


8x2  by1 


3x  3y 


or 


Therefore 


b"t 


at 


dZt 


at 


=  (A  +  2G)  V*  £ 


(15) 


(A  +  2G) 


VZ  £  -  V  2  •  V2  £ 


where 


-  / 


V  =  y  (A  +  2G)  /  p 

p 


This  is  the  wave  equation  and  shows  that  the  dilatation,  t 
(volumetric  strain)  is  propogated  through  the  medium  with  velocity, 
(X  +  2G)/p,  which  is  referred  to  as  compression  wave  or  P-wave. 


On  the  other  hand,  differentiating  both  sides  of  Equation  (13) 
with  respect  to  y,  and  of  Equation  (14)  with  respect  to  x, 
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a2 

<9u 

de 

du 

p 

(— )  - 

(  \  +  G  )  -  + 

GV2  (  —  )  . 

(16) 

„  2 

0  t 

ay 

3x3y 

ay 

and 

a2 

0  V 

3e 

dv 

p 

(  _  )  . 

(  A  +  G  )  -  + 

GV2  (  —  )  . 

(17) 

at2 

ax 

dxdy 

dx 

Subtracting 

Equation 

(16)  from  Equation 

(17)  yields, 

a2 

3u 

a  v 

0U 

)  =  G  V2  (  —  - 

av 

at2 

oy 

ax 

ay 

3x 

S  ince 

2w  =  - 

av 

8  u 

ax 

ay 

therefore , 


a2- 
a  u> 


=  G  V 


(18) 


at 


or 


a2- 

a  w 


v2  w  -  V  2  v2  U> 


(19) 


at 


where 


Equation  (19)  represents  the  equation  for  distortional  waves  and 
the  velocity  of  propogation  is  V  .  This  is  also  referred  to  as  the 
shear  wavt  or  S-Wave. 

We  now  consider  Rayleigh  wave.  Let  us  consider  a  plane  wave 
through  an  elastic  medium  with  a  plane  boundary,  and  let  u  and  v 
represent  the  displacements  in  the  x  and  y  directions,  respectively. 
Therefore , 


30  3\£ 

u  =  -  +  .  .  (20) 

dx  3y 


30 

84> 

3y 

dx 

where  0  and  t/>  are  two  potential  functions. 

The  dilation  e  can  be  defined  as 

3u  3v 

£  =  £  +  £  “  -  +  - 

x  y 

dx  dy 

or 

dZ0  3Z0 

£ - +  -  -  V2  0  . 

2  2 
dx  dy 

Similarly,  the  rotation  in  the  x-y  plane  can  be  given  by 


(21) 


(22) 


10 


3u  3v  d2 \J>  3  ij) 

2J>  =  -  -  -  =  -  +  -  -  V2V-  .  (23) 

3y  3x  3x  3y 

Substituting  Equations  (20)  and  (22)  into  Equation  (13)  yields, 

3  320  3  aV  8  2  8  z 

p  -  (  -  )  +  p  -  (  -  )  =  (A  +  2G)  —  (  V  0  )  +  G  (  V  V  ) 

3x  3t2  3y  3t2  3x  3y 

.  (24) 


Similarly,  substituting  Equations  (21)  and  (22)  into  Equation  (14) 
yields , 


3  a20  a 

^  -  (  -  )  -  p  — 

3y  3t2  3x 


a2^  a 

(  -  )  -  (  A  +  2G  )  - 

at2  ay 


a 

(  v20  )  -  g  —  (  V20  ) 

3x 

.  (25) 


Equations  (24)  and  (25)  will  be  satisfied  if  : 


(1) 


a20 


at' 


A  4-  2  G 

(  -  )  v2  0 

P 


V  2  •  V2 ■ 0 

P 


(2) 


a2v> 


at 


G 

-  .  v2  V  -  v  z-vz-\p 

p 


If  we  assume 
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and 


i  (wt 

F(y)  e 

fx) 

.  (26) 

N  i (wt  - 

C(y)  e 

fx) 

.  (27) 

where  F(y)  and  G(v)  are  functions  of  depths, 
f  -  2  -  7r  /  wavelength,  and 

i  -  y  - 1  ■ 

Solving  Equations  (26)  and  (27)  yields, 


-qv  i(wt  -  fx) 

0  -  A  •  e  •  e 


(28) 


„  -sy  i (wt  -  fx) 

xf>  =  B  •  e  J  •  e 


(29) 


where 


A  and  B  are  constants, 


(f2  -  w2  /  V  2) 

p 


(30) 


-/ 


(f2  -  w2  /  V  2) 


(31) 


Considering  boundary  conditions  of  ^  -  0  and  0  at  surface,  and 

combining  Equations  (20),  (21),  (28)  and  (29)  yields, 


16  (  1 


2  . 
U)  xp 


v  2f2 


)  (  i  - 


A  +  2G  w  w 

)  -  (  2  -  (  -  )  -  )2(  2 - )2 


2  -2 
V  f 


v2  f2 


v2f2 


(32) 
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2 


However,  V 

p 


=  (1/p)  (A  +  2G) 


and  V  2  - 

S 


G/p ,  therefore, 


2 

a  = 


s 


p 


G 


A  +  2G 


(35) 


In  Hooke's  law,  poisson's  ratio  is  defined  as 


1  A 

"  -  .  .  .  (36) 

2  (  A  +  G  ) 

From  Equation  (36), 

2  ■  i/  •  G 

*  =  .  .  (37) 

(  1  -  2v  ) 


Substituting  Equation  (37)  into  (35)  yields, 

(  1  -  2v  ) 

2 

Q  -  - 

(  2  -  2u  ) 

Again,  substituting  Equations  (33),  (34),  and  (  > 5)  into  (32) 
yields , 

V6  -  8V*  (16a2  -  24)  V2  -  16(  1  -  a  )  «=  0  .  (38) 


2 

This  equation  can  be  used  to  determine  the  proper  value  of  V  in  terms 

of  V  and  V  at  a  given  value  of  Poisson's  ratio.  Table  1  shows  some 
p  * 

values  of  V  -  V  /V  at  various  Poisson's  ratio, 
r  s 

For  displacement  of  the  Rayleigh  waves,  combining  and  solving 
Equations  (20),  (21),  (28),  and  (29)  yields, 
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u  =  Afi  (  -  e  ^  + 


and 


2qs 


2  c2 
s  +  f 


2qf 


-  s  y  .  i(wt  -  tx) 

e  )  e 


v  =  Aq  (-q/f  e  + - e  )  e 


2 

s  +  f 


TABLE  1.  VARIATION  OF  THE  RATIO  OF  THE  RAYLEIGH  TO 
THE  SHEAR  WAVE  FOR  VARIOUS  POISSON'S  RATIO 


V 

V  =  V  /V 

r  s 

0.25 

o.9i: 

0.29 

0.926 

0.33 

0.933 

0.40 

0.943 

o 

o 

0.955 

If  we  define  the  rate  of  attenuation  of  the  displacement  in 
directions , 


X(y) 


Y(y) 


- (q/f) (ty) 
e  + 


2(q/f)(s/f) 

(  -  ) 


2 

S  /f 


1 


-(s/f>(fy> 

e 


q  - (qf  t) (ty) 

-  e  + 

f 


2(q/f) 


2  /j-2  , 

S  /f  +  1 


■ (s/f ) (fy ) 

e 


(39) 


(40) 


x  and  y 


(41) 


(42) 
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For  example,  if  we  choose  ;<  poisson's  ratio,  v  -  0.25,  then 


X(y) 


Y(y) 


(-0.8*75  fy)  n  c  (-0.3933  fy) 

-e  +  0.5.73  e 

+0.8475  e‘-°  8475  fy)  -  1.4679  e1'03933  £y) 


(43) 

(44) 


Figure  3  shows  a  nondimens ional  plot  of  the  variation  of  the  amplitude 
of  vertical  and  horizontal  components  of  Raleigh  waves  with  depth  for 
various  poisson's  ratio. 


Based  on  Equations  (43)  and  (44)  and  Figure  3,  the  following 
observation  can  be  made. 


1.  The  magnitude  of  X  decrease  rapidly  with  increasing  value  of 

fy.  At  fv  «  1.21,  X  becomes  equal  to  zero;  so,  at  y 

2m/wavelength ,  thus  at  y  -  1.21/f  -  1 . 21 (wavelength)/2m  -  0.1V26 
(wavelength),  the  value  of  X  is  zero.  At  greater  depth,  X  becomes 
finite;  however  it  is  of  the  opposite  sign,  so  the  vibration  takes 
place  in  opposite  phase. 

2.  The  magnitude  of  Y  first  increases  with  fy,  reaches  a  maximum 

value  at  y  -  0.076  (wavelength)  (i.e,  fy  -  0.4775),  and  then 

decreases  with  depth. 


C.  EQUATION  OF  MOTION  IN  AN  ELASTIC  MEDIUM  WITH  GRAVITY  FORCE 

The  next  consideration  is  the  gravity  force  in  the  analysis.  On 
introducing,  g,  acceleration  of  gravity  in  the  y  direction,  into  the 
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Equations  (11)  and  (12),  we  obtain 


where  u  and  v  are  displacement  in  x  and  y  directions,  respectively,  and 
P  is  mass  of  density.  We  have  assumed  the  existence  of  a  uniform 
pressure  gradient  due  to  gravity  (i.e.,  the  initial  hydrostatic  pressure 
distribution).  Substituting  Hooke's  law  for  an  isotropic  medium 

a  =  Ac  +  2Gc 

X  X 

o  —Ac  +  2Gc 
y  y 

r  =  G  •  7 

xy  xy 

into  Equations  (45)  and  (46),  we  find 


(  A  +  G  ) 


+  GV  u  +  p  g 


(  A  +  G  ) 


+  GV  v  +  p  g 


By  differentiating  Equations  (47)  and  (48)  with  respect  to  x  and  y, 
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ft  spec  t  i  ve  1  v  ,  and  adding,  wu  obtain 


ail  civ 

t 

,’X  c'i  V 


O'  t  0 ‘  <  <3u  dv 

( x  + 1; )  (  —  +  —  )  +  gv7  (  —  ■<  —  ) 

dx2  dy2  dx  dy 


I 


d  u 


Pi\  (  -  “ 

dx2 


—  ) 
dy2 


V 

t/u 

—  -  ) 

and 

du 

;  .  ( - 

r  v 

X 

ti  v 

(  A  t  2G  )  V  i  +  pg  ST"  .  (49) 


( >n 

t.lu* 

other  hand 

re  s 

pt‘  <■  t 

to  v  and  >  , 

oh  t 

a  i  n 

bv  di  i  ferent  iat  ing  Equations  (47)  and  (48)  with 
reflectively,  and  subtracting  one  from  the  other  wo 


tiu 

(1  V 

dt 

)  -  (A  +  G)(  - 

a  c 

4 

O 

< 

bj 

i  £ 

i 

dv 

dy 

ax 

dxdy 

axay 

dy 

dx 

d 

+  pg  -  (  dv  -  dv  ) 

dx  dy 


or 


d  u 


P 


C,  V2w 


(50) 


Equations  (49)  differs  from  Equation  (15)  for  a  medium  without 
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gravity  force.  Equation  (49)  shows  the  existence  of  coupling  between 
longitudinal  and  transversal  waves.  This  coupling  is  due  to  the 
gravity  force.  The  physical  meaning  of  this  coupling  and  the 
significance  of  the  gravity  force  in  body  waves  should  be  investigated. 
A  comprehensive  study  will  be  carried  out  during  the  subsequent 
sections . 


20 


SECTION  III 


LABORATORY  INVESTIGATION 


A.  INTRODUCTION 

The  two-dimensional  governing  equations  of  motion  in  an 
elastic  medium  with  gravity  force  in  the  y  direction  were  given  in 
Equations  (45)  and  (46)  They  are 


da 


8x 


OT 


ov 


+  p  g 


dx 


32  u 


at 


(51) 


dr 


ax 


da 


xv  y 

— —  +  - — 


ay 


3v 


+  P  g 


ay 


a2v 


at 


(52) 


where  u  and  v  are  displacement  in  x  and  y  directions, 
respectively,  and  p  is  the  mass  density.  By  assuming  the  initial  state 
of  stress  is  hydrostatic  and  the  stress  -  strain  rel~*-ionship  for  an 
isotropic  medium  obeys  the  Hooke's  law,  we  obtain 


a2u  3«  3v 

p  -  -  (A  +  G)  -  +  G  V2u  +  pg  .  .  (53) 

at2  dx  dx 

a2v  3e  3v 

p  -  -  (A  +  G)  -  +  G  V2v  +  pg  .  .  (54) 

at2  ay  ay 
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Equations  (53)  and  (54)  can  be  simplified  by  an  additional 
assumption  to  gain  further  insight  of  the  roles  played  by  the  gravity 
terms,  namely  that  the  material  is  incompressible.  This  eliminates  the 
dilatational  wave.  Hooke ' s  law  can  then  be  written  for  an  incompressible 
material  as  follows: 


a  =  -p  +  2  •  G  •  e 

X  X 


(55) 


a  =  -p  +  2-G-f 

y  y 

t  =  G  •  7 

xv  xy 


(56) 

(57) 


where  p  is  the  increment  of  hydrostatic  pressure.  By  substituting, 
Equations  (55),  (56)  and  (57)  in  Equations  (51)  and  (52)  and  taking 

into  account  the  condition  of  incompressibility,  —  +  —  =  0, 

3x  3y 


we  find 


a2u  a 

P  -  =  G  V2u  -  -  (p  -  pgv)  .  (58) 

at2  3x 

a2v  a 

P  -  “  G  V2V  -  -  (p  -  pgv)  .  (59) 

3t2  3y 


If  one  anticipates  that  the  effect  of  gravity  is  important,  one  can 
derive  a  dimensionless  parameter,  which  represents  the  influence  of 
gravity,  from  Equation  (59).  This  requires 

3v 

G  •  V2  ■  v  -  pg  - 

3y 
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l .  e  .  , 


G-v  p- g-v  p- g-A 

-  -  or  - 

A2  A  G 


1 


where  A  is  the  wavelength  of  the  gravity  wave.  Therefore,  the 
dimensionless  gravity  parameter  r)=pgA/G  is  the  key  to  the  present 
investigation  of  gravity  effect.  After  defining  a  length  scale  J?=G/pg, 
the  gravity  parameter  rj  becomes  A/i .  In  general,  a  structure  is  most 
vulnerable  to  the  Rayleigh  wave  when  its  horizontal  length  scale  is 
close  to  half  of  the  wavelength.  For  instance,  £= 2,570  feet  for  a  soil 
with  G=3,000  lb/in2  and  pg=168  Ib/ft'3.  Suppose  that  the  gravity  effect 
is  important  when  r?>0 . 5 .  This  corresponds  to  A  >  1,285  feet  and  the 
horizontal  length  scale  of  the  structure  >  640  feet.  That  means  the 
structure  of  chat  length  is  vulnerable  to  the  Rayleigh  wave  enhanced  by 
gravity.  This  gravity  parameter,  r/=pgA/G,  will  be  varied  in  the 
laboratory  investigation  to  see  the  gravity  effects. 

Another  possible  gravity  effect  may  be  related  to  the  strength 
of  the  initial  blast  loading.  The  gravity  effect  is  significant  during 
the  initial  cratering  stage  after  explosion.  The  important  parameter 

3 

would  be  Q/(GA  )  ,  where  Q  is  the  impact  energy  transmitted  to  the 
soil . 


The  two  parameters  pgA/G  and  Q/GA  ,  identified  above,  will  be  key 
parameters  for  our  subsequent  laboratory  investigation . 
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B.  DESCRIPTION  OF  TESTING  FACILITY  AND  MODEL 

The  test  setup  consisted  of  a  drop  weight  mechanism,  a  model  box, 
and  electronic  equipments.  A  horizontal  steel  rod  that  could  be 
adjusted  to  any  desired  height  was  supported  by  a  rigid  steel  framework. 
A  smooth  1/8- inch  diameter  steel  rod  was  then  fixed  vertically  to  the 
horizontal  rod.  The  weight  was  free  to  slide  down  the  guide  when 
released,  and  the  motion  of  the  falling  weight  was  essentially 
equivalent  to  free  fall.  The  weight  was  cut  from  an  aluminium  stock  and 
machined  to  a  close  tolerance  to  allow  the  guide  to  pass  through  it 
smoothly . 

The  model  box  was  made  of  plexiglass,  18  inches  in  nominal  outside 
diameter  and  6  inches  high,  with  the  center  of  the  model  box  positioned 
directly  below  the  falling  weight.  The  material  used  in  the  model  tests 
was  commercially  available  gelatine.  The  gelatine  was  first  mixed  with 
hot  water  in  a  large  container,  then  was  carefully  poured  into  the 
plexiglas  model  box  up  to  5  inches  high.  The  model  was  placed  in  a  cool 
temperature  room  for  approximately  24  hours.  After  the  mixture  was 
hardened,  two  pressure  transducers  were  placed  on  the  surface  of  the 
medium:  one  at  the  center  and  the  other  one  at  a  point  6  inches  away 
from  the  center.  The  fluid  mixture  was  again  poured  into  the  model  box 
to  the  top  of  the  box,  and  the  model  was  placed  in  the  cool  temperature 
room  for  another  24  hours. 

The  pressure  transducers  used  in  the  tests  were  made  by  the 
scientists  at  the  General  Technology.  The  transducer  was  made  of 
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piezoelectric  ceramic  (see  Figure  4  )  .  The  vertical  surface  velocity 
was  measured  by  a  magnetic  velocity  transducer.  It  consists  of  a  small 
magnet  which  moves  freely  inside  a  coil.  The  motion  of  the  magnet  will 
induce  a  voltage  in  the  coil,  similar  to  an  electric  generator.  The 
transducer  was  carefully  calibrated  against  known  velocities  for 
quantitative  measurements. 

Two-channel,  type  MARK220,  special  purpose  recorder  manufactured  by 
the  Brush  Instruments  was  used  to  record  the  pressure  induced  in  the 
model  by  the  falling  weight.  All  data  before  and  after  impact  loading 
was  recorded  and  printed  out  on  a  hard  copy.  A  photograph  of  the  drop 
weight  system,  model  box  filled  with  the  mixture  (gelatine),  embedded 
pressure  transducers  and  two-channel  recorders  is  shown  in  Figure  5. 

Seismic  methods,  including  the  crosshole,  downhole,  surface 
refraction,  surface  reflection,  steady-state  Rayleigh-wave  method  and 
the  spectral-analysis-of-surface-waves(SASW)  method,  are  often  used  to 
determine  the  elastic  material  properties  (i.e.,  shear  modulus,  etc.). 
However,  in  this  experimental  study,  the  shear  modulus  of  the  model 
sample  was  determined  by  a  different  method  that  was  less  time  consuming 
and  complicated  than  the  seismic  methods.  Rather,  wc  used  a  method 
simple,  yet  accurate  enough  to  measure  the  shear  modulus  of  the  model. 
A  cube  sample  with  dimensions  of  1  inch  x  1  inch  x  1  inch  was  made  of 
gelatine.  The  bottom  of  the  sample  was  fixed,  and  an  incremental 
shear  force  was  applied  on  the  top  surface  of  the  sample  by  pulling  it 
horizontally  through  a  frictionless  pulley.  The  horizontal  small 
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displacement  was  measured  through  a  microscope,  and  corresponding  strain 


was  calculated.  A  photograph  of  the  shear  loading  apparatus  and 
microscope  is  shown  in  Figure  6. 

C .  LABORATORY  PROGRAM 

This  investigation  assumed  that  the  vertical  movement  of  the  medium 

surface  due  to  impact  loading  of  the  drop  weight  would  be  important. 

3 

The  impact  parameter  is  defined  as  Q/(GA  )  where  Q  -  WH  is  the  impact 
energy  transmitted  to  the  soil,  W  is  the  drop  weight,  H  is  the  drop 
height,  A  is  the  wavelength  of  the  Rayleigh  wave  and  G  is  the  shear 
modulus  of  the  soild.  Two  different  masses  (0.2  pound  and  0.4  pound) 
and  three  different  diameters  (1  inch.,  2  inches.,  and  3  inches.)  of  the 
drop  weight  (see  Figure  7),  and  six  different  drop  heights  (1  inch,  2 
inches,  3  inches,  4  inches,  5  inches,  and  6  inches.)  were  chosen  as 
parameters  for  the  experimental  model  study.  Another  important 
parameter  to  match  between  the  laboratory  model  and  the  prototype  is  the 
gravity  parameter,  pgX/G.  This  is  a  measure  of  the  gravity  potential 
energy  in  relation  to  the  strain  energy  of  the  solid  medium.  To  study 
the  significant  gravity  effects,  it  is  essential  that  the  gravity 
parameter  pgA/G  be  kept  at  about  order  one.  To  achieve  this,  the  shear 
modulus  of  the  model  must  be  small  because  of  the  limited  size  of  the 
model.  As  a  result,  two  shear  modulii  of  0.2  and  1.4  psi  were  chosen 
for  the  laboratory  model  study.  Table  2  shows  the  parameters  considered 
in  the  model  tests. 
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Photograph  of  Shear  Loading  Apparatus 


TABLE  2.  PARAMETERS  CONSIDERED  IN  THE  MODEL  TESTS 


Parameter 

Shear  Modulus  of 

Mode  1 

0.2  psi  and 

1.4 

psi 

Drop  Weight 

0.2  lb .  and 

0.4 

lb. 

Diameter  of  Drop 

Weight 

1  in. ,  2  in. , 

and 

3  in . 

Drop  Height 

1  in . ,  2  in . ,  3 

in.  , 

,  4  in . , 

5  in . ,  and  6  in . 
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SECTION  IV 


LABORATORY  RESULT 

A.  THE  INFLUENCE  OF  SHEAR  MODULUS  OF  SOLID 

Seventy  two  model  tests  were  performed,  and  the  pressure  and 
velocity  responses  of  the  medium  under  different  loading  conditions  were 
obtained.  Figures  8  and  9  show  the  typical  responses  of  the  pressure 
and  velocity  at  1  inch  drop  height  with  a  3  inch  diameter  and  0.4  pound 
of  the  drop  weight  for  G  =  0.2  psi  and  G  —  1.4  ps: , respectively .  Based 
on  the  data  obtained,  frequency,  f  and  period,  T,  were  calculated, 
and,  accordingly,  the  wave  speed,  c,  and  wavelength,  A,  were  calculated 
by  the  following  equations  : 

w  =  2 7T  (-) 

T 

c  -  L/T 
A  =  c-T 

where  , 

w  is  the  frequency  of  the  harmonic  in  rad/sec, 

T  is  the  period  in  second, 

L  is  the  travel  distance  in  inch, 
t  is  travel  time  in  second, 
c  is  the  wave  speed  in  in/sec,  and 
A  is  the  wavelength  in  inch. 

The  period,  T  was  determined  by  finding  the  length  at  which  the  waves 
makes  a  complete  cycle .  A  few  of  these  cycles  make  up  a  group  of  waves , 
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therefore,  the  period  found  was  the  average  of  all  these  values.  After 
the  wavelength  has  been  determined,  two  dimensionless  parameters  were 
calculated : 

gravity  parameter  -  p gA/G 

3 

impact  parameter  -  W-H/(GA  ) 

where , 

p  is  density  of  the  medium  in  lb-sec2/in<' , 

2 

g  is  the  acceleration  of  the  gravity,  1  g  -  386  in/sec  , 

G  is  the  shear  modulus  of  the  medium  in  psi, 

W  is  the  drop  weight  in  pound,  and 
H  is  the  height  of  the  drop  weight  in  inch. 

Tables  3  thru  14  summarize  the  test  results.  As  shown  in  Tables  6 
and  12,  the  models  failed  after  2-inch  drop  tests.  Cracks  were 
developed  in  the  model  from  high  stress  concentration  caused  by  heavier 
weight  with  small  diameter  of  the  drop  weight.  Also,  as  shown  in  Tables 
9  and  11,  insufficient  data  were  obtained  during  the  model  tests  due  to 
electrical  disconnection  in  the  data  acquisition  system. 

Based  on  the  tables  obtained,  10  figures  were  drawn.  Figures  10 
and  11  represent  a  relationship  between  wave  speed  and  drop  height  of 
the  weight  for  models  with  G-0.2  psi  and  G-1.4  psi,  respectively.  It 
can  be  observed  that  wave  speed  of  each  model  remains  almost  constant  at 
different  drop  heights,  although  with  slight  variation.  Therefore,  the 
wave  speed  is  not  a  function  of  the  drop  height.  However,  the  wave  speed 
is  slightly  dependent  on  the  diameter  of  drop  weight.  The  3- inch 
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TABLE  3.  TEST  RESULTS  FOR  MODEL  1 


Shear  Modulus,  G  =  0.2  psi 
Drop  Weight,  W  =  0.2  lb. 
Diameter  of  Weight,  D  =  1.0  in. 


Drop  Height  (H) 

H  =  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

Model  (p) 

0.906x10"'' 

( lb-  seedin'1) 

0. 906x10"'' 
(lb-seedin'1) 

0.906x10"'’ 

(lb-seedin'”) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

36.3  in/sec 

39.0  in/sec 

38.0  in/sec 

Wavelength  (A) 

5.10  in 

5.50  in 

5.41  in 

Gravity  Parameter 
(pgVG) 

0.893 

0.963 

0.947 

Impact  Parameter 
(WH/G a  ) 

0.0075 

0.012 

0.0190 

Maximum  Pressure 

23.9  psi 

40.2  psi 

51.1  psi 

Maximum  Velocity 

2 . 50in/sec 

3 . 50in/sec 

3 . 90in/ sec 

Drop  Height  (H) 

H  =  4  in. 

H  =  5  in. 

H  =  6  in . 

Density  of 

o.  906x10“ 

0.906x10"'' 

0. 906x10“' 

Model  (p) 

2  4 

(lb-sec  /in  ) 

2  4 

(lb-sec  /in  ) 

(lb-seedin'1) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

37.3  in/sec 

38.5  in/sec 

40.6  in/sec 

Wavelength  (A) 

5.10  in 

5.45  in 

5.46  in 

Gravity  Parameter 
(PgVG) 

0.893 

0.954 

0.956 

Impact  Parameter 
(WH/GA3) 

0.0301 

0.0309 

0.0370 

Maximum  Pressure 

54.4  psi 

87.0  psi 

95.7  psi 

Maximum  Velocity 

3 . 80in/sec 

4 . OOin/sec 

4 . 40in/sec 
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TABLE  4.  TEST  RESULTS  FOR  MODEL  2 


Shear  Modulus,  G  -  0.2  psi 
Drop  Weight,  W  -  0.2  lb. 
Diameter  of  Weight,  D  —  2.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  =  3  in. 

Density  of 

Model  (/>) 

0 . 906x10" * 
(lb-sec2/inA) 

0 . 906xl0A 
(lb-sec2/inA) 

0 . 906x10"A 
(lb-sec2/inA) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

40.1  in/sec 

43.0  in/sec 

40.4  in/sec 

Wavelength  (A) 

4.20  in 

4.25  in 

4.49  in 

Gravity  Parameter 
(pgA/G) 

0.7350 

0.7439 

0.7860 

Impact  Parameter 
(WH/GAJ) 

0.0135 

0.0261 

0.0220 

Maximum  Pressure 

17.4  psi 

23.4  psi 

30.5  psi 

Maximum  Velocity 

3 . 00in/sct- 

3 . 90in/sec 

4 . 40in/sec 

Drop  Height  (H) 

H  =  4  in. 

H  -  5  in. 

H  -  6  in. 

Density  of 

0 . 906x10’A 

0. 906xl0"A 

0 . 906xl0"A 

Model  ( p ) 

(lb-sec2/inA) 

(lb- sec2/inA) 

(lb-sec2/inA) 

Specific  Gravity 

0.9694 

0.9694 

0.9694 

of  Model 

Wave  Speed  (c) 

40.3  in/sec 

42.4  in/sec 

43.0  in/sec 

Wavelength  (A) 

4.48  in 

4.70  in 

4.60  in 

Gravity  Parameter 

0.7840 

0.8228 

0.805 

(P  gVG) 

Impact  Parameter 

0.0445 

0.0482 

0.0616 

(WH/GA  ) 

Maximum  Pressure 

39.2  psi 

51.1  psi 

54.4  psi 

Maximum  Velocity 

4 . 95in/sec 

5 . lOin/sec 

5 . 35in/sec 
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TABLE  5.  TEST  RESULTS  FOR  MODEL  3 


Shear  Modulus,  G  -  0.2  psi 
Drop  Weight,  W  -  0.2  lb. 
Diameter  of  Weight,  D  -  3.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

Model  (p) 

0 . 906x10" 4 
(lb-sec2/in*) 

0.906x10'* 

(lb-s‘='c2/in*) 

0.906x10'* 

2  4 

(lb-sec  /in  ) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

45.2  in/sec 

45.3  in/sec 

45.7  in/sec 

Wavelength  (A) 

4.38  in 

4.39  in 

4.37  in 

Gravity  Parameter 
(pgA/G) 

0.7670 

0.7680 

0.7650 

Impact  Parameter 
(WH/GA  j 

0.0119 

0.0236 

0.0359 

Maximum  Pressure 

7.80  psi 

9.50  psi 

11.9  psi 

Maximum  Velocity 

2 . 80in/sec 

3 . 80in/sec 

4 . lOin/sec 

Drop  Height  (H) 

H  =  4  in. 

H  -  5  in. 

H  =  6  in. 

Density  of 

o.906xio'* 

0.906x10'* 

0.906x10'* 

Model  (p) 

( lb-sec2/in*) 

(lb-secZ/in*) 

(lb-secZ/in*) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

45.6  in/sec 

45.3  in/sec 

45.6  in/sec 

Wavelength  (A) 

4.35  in 

4.37  in 

4.42  in 

Gravity  Parameter 
(PgA/G) 

0.761 

0.7650 

0.774 

Impact  Parameter 
(WH/GA3) 

0,0486 

0.0599 

0.069 

Maximum  Pressure 

14.6  psi 

15.2  psi 

19.0  psi 

Maximum  Velocity 

4 . 70in/sec 

4 . 90in/sec 

5 . 20in/sec 
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TABLE  6.  TEST  RESULTS  FOR  MODEL  4 

Shear  Modulus,  G  -  0.2  psi 
Drop  Weight,  W  -  0.4  lb. 
Diameter  of  Weight,  D  -  1.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

0 . 906xl0’4 

0. 906x10"* 

0. 906xl0"A 

Model  (p) 

(lb-sec2/in**) 

(lb-secVin*) 

(lb-sec2/in<’) 

Specific  Gravity 

0.9694 

0.9694 

0.9694 

of  Model 

Wave  Speed  (c) 

42.0  in/sec 

48.7  in/sec 

* 

Wavelength  (A) 

6.56  in 

7.51  in 

* 

Gravity  Parameter 

1.185 

1.314 

* 

(pgA/G) 

Impact -Parameter 

0.0032 

0.0047 

* 

(WH/GA  ) 

Maximum  Pressure 

29.4  psi 

52.2  psi 

* 

Maximum  Velocity 

2 . 70in/sec 

4 . 33in/sec 

* 

Drop  Height  (H) 


H  -  4  in. 


H  -  5  in. 


H  *  6  in . 


Density  of 

Model  ( p ) 

Specific  Gravity 


0.906x10 


-4 


0.906x10 


0.906x10 


(lb-sec2/'in*)  (lb-secVin*)  (lb-se^/in^) 


0.9694 


0.9694 


0.9694 


of  Model 

Wave  Speed  (c) 

* 

* 

* 

Wavelength  (A) 

* 

* 

* 

Gravity  Parameter 

* 

* 

* 

(PgA/G) 

Impact  Parameter 

* 

* 

* 

(WH/GA  ) 

Maximum  Pressure 

* 

* 

* 

Maximum  Velocity 

* 

* 

* 

*  :  Model  Failure 
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TABLE  7.  TEST  RESULTS  FOR  MODEL  5 


Shear  Modulus,  G  -  0.2  psi 
Drop  Weight,  W  -  0.4  lb. 
Diameter  of  Weight,  D  -  2.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

’  -  ,  *  -  -  -c 

Model  ( p ) 

o  oncvid'* 

( lb  -  sec2/in“) 

0. 90^x10 "A 
(lb-sec2/in“) 

0.906x10"“ 

( lb- sec2/in“) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

46 . 2  in/sec 

46.1  in/sec 

46.9  in/sec 

Wavelength  (A) 

5.92  in 

5.92  in 

5.91  in 

Gravity  Parameter 

(p  gVG) 

1.036 

1.0340 

1.034 

Impact  Parameter 
(WH/GA-3) 

0.0096 

0.0194 

0.0290 

Maximum  Pressure 

25.6  psi 

45.7  psi 

59.1  psi 

Maximum  Velocity 

3 . 50in/sec 

4 . 80in/sec 

5 . 40in/sec 

Drop  Height  (H) 

H  =  4  in. 

H  -  5  in. 

H  -  6  in. 

Density  of 

0 . 906x10’* 

0.906x10"“ 

0.906x10"“ 

Model  ( p ) 

(lb-sec2/in“) 

(lb-sec2/in“) 

(lb-sec2/in“) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

46 . 2  in/sec 

46.6  in/sec 

46.6  in/sec 

Wavelength  (A) 

5.91  in 

5.87  in 

5.91  in 

Gravity  Parameter 
(pgA/G) 

1.033 

1.027 

1.034 

Impact  Parameter 
(WH/GA3) 

0.0390 

0.0494 

0.0582 

Maximum  Pressure 

80.5  psi 

135.9  psi 

163.2  psi 

Maximum  Velocity 

5 . 90in/sec 

6 . 80in/sec 

7 . 70in/sec 
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TABLE  8.  TEST  RESULTS  FOR  MODEL  6 


Shear  Modulus,  G  -  0 . 2  psi 
Drop  Weight,  W  -  0.4  lb. 
Diameter  of  Weight,  D  -  3.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

c;  i  £V  O  f 

Model  (p) 

0 . 906xl0’4 

(lb-  secVin*) 

0.906x1  o' 4 

(lb-secVin*) 

0 . 906x10*' 
(lb-sec^/in**) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

48.2  in/sec 

48 . 1  in/sec 

48.4  in/sec 

Wavelength  (A) 

4.44  in 

4.44  in 

4.43  in 

Gravity  Parameter 
(pgX/G) 

0.7772 

0.777 

0.7754 

Impact  Parameter 
(WH/GA  j 

0.0456 

0.0460 

0.1380 

Maximum  Pressure 

14.7  psi 

23.0  psi 

27.2  psi 

Maximum  Velocity 

3 . 90in/sec 

5 . 60in/sec 

5 . 90in/sec 

Drop  Height  (H) 

H  =  4  in. 

H  -  5  in. 

H  -  6  in. 

Density  of 

o^oexio"1 

o-gooxio"1 

0 . 906xl0"4 

Model  ( p ) 

(lb-se^/in4) 

(lb-sec2/inA) 

(lb-sec2/in*) 

Specific  Gravity 
of  Model 

0.9694 

0.9694 

0.9694 

Wave  Speed  (c) 

48.2  in/sec 

48.3  in/sec 

48.2  in/sec 

Wavelength  (A) 

4.44  in 

4.44  in 

4.44  in 

Gravity  Parameter 
(pgA/G) 

0.777 

0.7772 

0.777 

Impact  Parameter 
(WH/GA3) 

0.0916 

0.2284 

0.1370 

Maximum  Pressure 

34.8  psi 

38.1  psi 

48.9  psi 

Maximum  Velocity 

7 . lOin/sec 

7 . 70in/sec 

8 . 20in/sec 

TABLE  9.  TEST  RESULTS  FOR  MODEL  7 


Snear  Modulus,  G  -  1.4  psi 
Drop  Weight,  W  -  0.2  lb. 
Diameter  of  Weight,  D  «  1.0  in. 


Drop-Height  (H) 

H  =  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

0 . 895x10" 4 

0 . 895x10' A 

0 . 895x10" ^ 

2  .  4 . 

Model  (p) 

viu  - /  j.n  ) 

(lb  sec  /in  N 

f  -  cor  /itl  i 

Specific  Gravity 

0.9577 

0.9577 

0.9577 

of  Model 

Wave  Speed  (c) 

105 . 5in/sec 

106 . 4in/sec 

110 . lin/sec 

Wavelength  (A) 

5.52  in 

5.60  in 

5.94  in 

Gravity  Parameter 

0.135 

0.138 

0.147 

(pgVC) 

Impact -Parameter 

0.0008 

0.0016 

0.002 

(WH/GA  ) 

Maximum  Pressure 

60.9  psi 

80.1  psi 

100.1  psi 

Maximum  Velocity 

1 . 26in/sec 

1 . 90in/sec 

3 . OOin/sec 

Drop  Height  (H)  H  -  4  in.  H  =  5  in.  H  ■=  6  in. 


Density  of 

o.sgsxio''1 

0 . 895xl0"4 

0.895x10'* 

Model  (p) 

2  4 

(lb-sec  /in  ) 

2  4 

(lb-sec  /in  ) 

(lb-sec2/:.!-/) 

Specific  Gravity 

0.9577 

0.9577 

0.9577 

of  Model 

Wave  Speed  (c) 

113 . lin/sec 

116 . 7in/sec 

** 

Wavelength  (A) 

6.50  in 

5.92  in 

** 

Gravity  Parameter 

0.161 

0.146 

** 

(PgVC) 

Impact  Parameter 

0.0021 

0.0034 

** 

(WH/GA  ) 

Maximum  Pressure 

135.9  psi 

160.1  psi 

** 

Maximum  Velocity 

3 . 60in/sec 

3 . 80in/sec 

** 

**  :  Insufficient 

Data 
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TABLE  10.  TEST  RESULTS  FOR  MODEL  8 


Shear  Modulus,  G  =  1.4  psi 
Drop  Weight,  W  =  0.2  lb. 
Diameter  of  Weight,  D  =  2.0  in. 


Drop  Height  (H) 

H  =  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

0.895x10"'’ 

0.895x10"'’ 

O^SkIO"’ 

1  / 

(  1  W  .  ,-r.^2  /i 

(lv  -  sec2/ in'") 

( lh- sec.  /i n  ) 

Specific  Gravity 

0.9577 

0.9577 

0.9577 

of  Mod<=l 

Wave  Speed  (c j 

117 . 9in/sec 

116 . 8in/sec 

118 . 9in/sec 

Wavelength  (A) 

5.42  in 

5.^3  in 

5.42  in 

Gravity  Parameter 

0.1339 

0.1342 

0.1339 

(pgA/G) 

Impact Parameter 

0.0010 

0.0018 

0.0028 

(WH/GA  ) 

Maximum  Pressure 

ZZ .  5  psi 

78.9  psi 

90 .30  psi 

Maximum  Velocity 

2 . 50in/sec 

3 . 40in/sec 

3 . 50in/sec 

Drop  Height  (H) 

H  =  4  in. 

H  -  5  in. 

H  =  6  in. 

Density  of 

0.895x10'* 

0.895x10"'’ 

0.895x10'* 

Model  ( p ) 

(lb-sec2/in*) 

(lb-secVin'*) 

(lb-sec2/in'‘) 

Specific  Gravity 

0 . 9577 

0.9577 

0.9577 

of  Model 

Wave  Speed  (c) 

118 . 9in/sec 

118 . 5in/sec 

118 . 5in/sec 

Wavelength  (A) 

5.42  in 

5.35  in 

5.30  in 

Gravity  Parameter 

0.1339 

0.1320 

0.131 

(P gA/G) 

Impact  Parameter 

0.0036 

0.0047 

0.0058 

(WH/GA3) 

Maximum  Pressure 

104.5  psi 

123.1  psi 

155.4  psi 

Maximum  Velocity 

3 . 85in/sec 

4 . 30in/sec 

4 . 75in/sec 

42 


TABLE  11.  TEST  RESULTS  FOR  MODEL  9 


Shear  Modulus,  G  -  1.4  psi 
Drop  Weight,  W  —  0.2  lb. 
Diameter  of  Weight,  D  -  3.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

Mod^1  f  " ' 

0.895x10''' 

( lb-sec2/in<‘) 

0.895x10** 

2  4 

(lb -sec  /in  ) 

0.895x10"'' 
(lb- seedin'1) 

Specific  Gravity 
of  Model 

0.9577 

0.9577 

0.9577 

Wave  Speed  (c) 

122 ,4in/sec 

** 

123 . 3in/sec 

Wavelength  (A) 

4.53  in 

** 

5.90  in 

Gravity  Parameter 
(pgA/G) 

0.1006 

** 

0.1313 

Impact  Parameter 
(WH/GA  ; 

0.0016 

** 

0.0020 

Maximum  Pressure 

25.0  psi 

** 

73.9  psi 

Maximum  Velocity 

2 . 80in/sec 

** 

5 . 17in/sec 

Drop  Height  (H) 

H  -  4  in. 

H  -  5  ir. 

H  -  6  in. 

Density  of 

0 . 895x10"'' 

0. 895x10"'’ 

0 . 895x* 

Model  ( p ) 

(lb-sec2/in<') 

(lb -seedin'1) 

(lb-sec2/in'1) 

Specific  Gravity 
of  Model 

0.9577 

0.9577 

0.9577 

Wave  Speed  (c) 

** 

123 . 3in/sec 

** 

Wavelength  (A) 

** 

5.90  in 

** 

Gravity  Parameter 
(pgA/G) 

** 

0.1095 

** 

Impact  Parameter 
(WH/GA3) 

** 

0.0060 

** 

Maximum  Pressure 

** 

78.3  psi 

** 

Maximum  Velocity 

** 

6.61  in/sec 

** 

** 


Insufficient  Data 
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TABLE  12.  TEST  RESULTS  FOR  MODEL  10 


Shear  Modulus,  G  -  1.4  psi 
Drop  Weight,  W  -  0.4  lb. 
Diameter  of  Weight,  D  -  1.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

Model  (p) 

0.895x10"* 

2  4 

(lb- sec  /in  ) 

0.895x10"* 

(lb-secZ/in*") 

0.895x10'* 

(lb-secZ/in*) 

Specific  Gravity 
of  Model 

0.9577 

0.9577 

0.9577 

Wave  Speed  (c) 

192 . 5in/sec 

181 . 7in/sec 

* 

Wavelength  (A) 

7.89  in 

6.80  in 

* 

Gravity  Parameter 
(pgA/G) 

0.195 

0.168 

* 

Impact  Parameter 
(WH/GA  ; 

0.0003 

0.0009 

* 

Maximum  Pressure 

71.8  psi 

122.3  psi 

* 

Maximum  Velocity 

2 . 16in/sec 

3 . 25in/sec 

* 

Drop  Height  (H)  H  =  4  in.  H  =  5  in.  H  -  6  in. 

Density  of  0.895x10"*  0.895x10"*  0.895x10"* 

Model  (p)  (lb-secZ/in*)  (lb-secZ/in*)  (lb-sec2/in* 

Specific  Gravity  0.9577  0.9577  0.9577 

of  Model 


Wave  Speed  (c) 

★ 

* 

* 

Wavelength  (A) 

* 

* 

* 

Gravity  Parameter 

* 

* 

* 

(pgVG) 

Impact  Parameter 

* 

* 

* 

(WH/GA  ) 

Maximum  Pressure 

★ 

★ 

* 

Maximum  Velocity 

* 

★ 

* 

*  -  Model  Failure 
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TAL'Li:  13.  TEST  RESULTS  FOR  MODEL  11 


Shear  Modulus,  G  -  1.4  psi 
Drop  Weight,  W  =  0.4  lb. 
Diameter  of  Weight,  D  =  2.0  in. 


Drop  Height  (H) 

H  =  1  in. 

H  -  2  in. 

H  =  3  in. 

Density  of 

Model  (p) 

0 . 895x10'" 

2  4 

(lb-sec  /in  ) 

0.895x10'" 

(lb-sec2/in") 

0.895x10'" 

7  *4 

(lb-sec  /in  ) 

Specific  Gravity 
of  Model 

0.9577 

0.957/ 

0.9577 

Wave  Speed  (c) 

114 . 8in/sec 

114 . 8in/sec 

114 . 8in/sec 

Wave  length  (A) 

4.44  in 

4.85  in 

4.80  in 

Gravity  Parameter 
(pgA/G) 

0.109 

0.1198 

0.1186 

Impact  Parameter 
(WH/GA; 

0.0032 

0.0050 

0.0078 

Maximum  Pressure 

67.5  psi 

95.7  psi 

110.3  psi 

Maximum  Velocity 

4 . lin/sec 

4 . 90in/sec 

5 . 30in/sec 

Drop  Height  (H) 

H  ■=  4  in . 

H  -  5  in. 

H  =  6  in. 

Density  of 

0.895x10'" 

0 . 895x10'* 

0.895x10'" 

Model  (p) 

2  4 

(lb-sec  /in  ) 

(lb-sec2/in") 

(lb -sec  /in") 

Specific  Gravity 
of  Model 

0.9577 

0.9577 

0.9577 

Wave  Speed  (c) 

116 . 8in/sec 

116 . 9in/sec 

117 . Oin/sec 

Wavelength  (A) 

4.46  in 

4.54  in 

5.00  in 

Gravity  Parameter 
(PgA/G) 

0.110 

0.112 

0.124 

Impact  Parameter 
(WH/GA3) 

0.0129 

0.0150 

0.0137 

Maximum  Pressure 

150.5  psi 

184.9  psi 

230.0  psi 

Maximum  Velocity 

5 . 9in/sec 

6 . 6in/sec 

7 . 2in/sec 

TABLE  14.  TEST  RESULTS  FOR  MODEL  12 


Shear  Modulus,  G  -  1.4  psi 
Drop  Weight:,  W  —  0.4  lb. 
Diameter  of  Weight,  D  -  3.0  in. 


Drop  Height  (H) 

H  -  1  in. 

H  -  2  in. 

H  -  3  in. 

Density  of 

Model  (p) 

0.895x10'* 

2  4 

(lb-sec  /in  ) 

0.895x10'* 

(lb-secVin*1) 

0.895x10'* 

(lb-sec2/in*) 

Specific  Gravity 
of  Model 

0.9577 

0.9577 

0.9577 

Wave  Speed  (c) 

120 . lin/sec 

120 . lin/sec 

121 . lin/sec 

Wavelength  (A) 

4.99  in 

4.99  in 

4.95  in 

Gravity  Parameter 

(Pg>/C) 

0.1233 

0.123 

0.1223 

Impaet,r^ra.r.eter 
( WH/GA  ; 

0.0044 

0.0046 

0.014 

Maximum  Pressure 

46.8  psi 

60.2  psi 

82.4  psi 

Maximum  Velocity 

4 . 6in/sec 

5 . 65in/sec 

7 . 2in/sec 

Drop  Height  (H) 

H  -  4  in. 

H  -  5  in. 

H  «=  6  in. 
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Figure  1]  .  Wavelength  vs.  Drop  Height  for  Models  with  G=1.4  psi 


diameter  dt_op  weight  induces  the  highest  wave  speed  while  the  1  -  inch 
diameter  drop  weight  induces  the  lowest  wave  speed,  and  may,  thus 
indicate  the  influence  of  cratering  sizes  on  Rayleigh's  wave  speed. 
However,  the  size  of  the  cratering  is  coupled  with  impact  and  explosion 
energy,  thereby,  impact  loadings  and  explosion  energies  are  important 
factors  in  practice.  As  expressed  earlier,  a  higher  wave  speed  was 
obtained  when  a  higher  shear  modulus  (1.4  psi)  of  the  medium  was  used. 
It  confirms  the  influence  of  shear  modulus  on  propogating  wave  speed. 

Figures  12  and  13  show  a  relationship  between  wavelength  and  drop 
height  of  the  drop  weight  for  models  with  G  =  0.2  psi  and  G  -  1.4  psi, 
respectively.  First  observation  from  the  figure  again,  shows  the 
independency  of  the  wavelength  from  the  drop  height.  Second  observation 
is  that  the  magnitude  of  the  wavelength  are  between  4  to  6  inches  for 
all  tests.  However,  it  is  very  difficult  to  define  the  parameters  that 
influence  the  differences  in  wavelength  magnitude.  It  is  certainly  due 
to  period,  T,  and  wave  speed  which  is,  in  turn,  related  to  the  diameter 
of  the  drop  weight. 

Figures  14  and  15  represent  a  relationship  between  two 
dimensionless  parameters,  pgA/G  and  pgD/G,  for  models  with  G  -  0.2  psi 
and  G  -  1.4  psi.  The  wavelength  is  dependent  on  the  diameter  of  the 
drop  weight;  however,  this  dependence  on  the  diameter  is  very  slight 
when  W  -  0.2  pound.  It  also  shows  that  the  dependence  cn  the  drop 
height  is  insignificant. 

The  maximum  pressure  recorded  by  the  pressure  transducer  at  the 
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Figure  ]_2.  Wave  Speed  vs.  Drop  Height  for  Models  with  G=0.2  psi 
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Figure  15.  A  Relationship  of  Dimensionless  Parameters  for  Models 


center  represents  the  degree  of  impact  as  well  as  the  magnitude  of 
deformation  at  that  particular  location.  Another  indication  of  the 
magnitude  of  deformation  of  the  medium  caused  by  impact  loading  is 
represented  by  the  maximum  vertical  velocity  recorded  by  the  off-center 
velocity  transducer  at  the  surface. 

Figures  16  and  17  show  the  relationship  between  the  maximum 
pressure,  P  and  the  impact  energy,  WH .  The  results  show  that  the 
maximum  pressure  has  a  linear  relationship  with  the  impact  energy,  WH , 
and,  as  expected,  the  smaller  diameters  induce  higher  maximum  pressures 
at  a  given  energy  due  to  the  smaller  contact  area. 

Figures  18  and  19  show  the  relationship  between  the  maximum 
vertical  velocity  of  the  medium  on  the  surface  and  the  impact  energy, 
WH ,  for  both  cases. G  =  0.2  psi  and  1.4  psi.  As  shown  in  Figures  10  and 
11  for  the  wavespeed,  the  vertical  velocity  is  coupled  with  the  diameter 
of  the  drop  weight.  The  3 -inch  diameter  drop  weight  induces  the  highest 
vertical  velocity  while  1-inch  diameter  drop  weight  induces  the  lowest 
vertical  velocity  for  all  cases.  The  results  also  show  the  linear 
dependence  on  WH  for  most  cases.  The  magnitudes  of  the  velocity  vary 
from  1  to  9  in/s. 
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Figure  16.  Maximum  Pressure  vs.  Impact  Energy  for  Models  with  G=0.2  psi 
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Figure  18.  Maximum  Velocity  vs.  Impact  Energy  for  Models  with  G-0.2  psi 


56 


Figure  29.  Maximum  Velocity  vs.  Impact  Energy  for  Models  with  G=l. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  following  conclusions  were  drawn  from  the  analysis  of  the 
equation  of  motion  in  an  elastic  medium  with  gravity  force  and  the 
laboratory  investigation  performed. 

1.  Wave  speed,  c,  is  constant  and  independent  of  the  drop  height, 
H,  at  a  given  diameter  of  the  drop  weight.  However,  wave  speed 
increases  slightly  when  the  diameter  of  the  drop  weight,  W,  increases. 

2.  The  wavelength  is  virtually  independent  of  the  drop  height  at 
a  given  diameter  of  the  drop  weigh*".  The  wavelength  is  decreased  while- 
diameter  of  the  weight  increases.  This  is  shown  contradictory  to  our 
intuition  that  the  wavelength  would  be  larger  due  to  larger  contact 
areas.  Further  investigation  is  needed  for  clearer  understanding  on  this 
phenomenon . 

3.  The  governing  parameter  in  this  study  is  the  gravity  parameter, 
p gA/G.  The  gravity  effect  is  prominent  when  the  gravity  parameter  is  of 
order  of  1  or  larger.  It  is  a  function  of  the  intensity  of  ground 
movement  which  is  related  to  the  intensity  of  vertical  displacement  due 
to  impact  loadings.  When  gravity  effect  is  important,  the  scaling 
length  for  ground  movement  is  -?-=G/p g. 
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4.  Gravity  parameter.  pgA/G,  varies  from  0.72  to  1.00  for  a  model 
with  shear  modulus,  G  =  0.2  psi  and  from  0.12  to  0.14  for  a  model  with 
G  =  1.4  psi.  It  indicates  that  the  use  of  a  model  with  large  shear 
modulus  is  not  adequate  to  investigate  the  gravity  effects  since  the 
model  is  too  stiff,  thereby ,  small  surface  movement  is  induced. 

5.  The  governing  parameter  pgA/G  is  uniquely  related  to  another 
parameter  pgD/G  under  impact  loading.  Therefore,  there  would  be  a 
relationship  between  the  wave  length,  A  and  the  diameter  of  the 
drop  weight,  D.  Further  investigation  is  needed  to  study  these 
relationships . 

6.  The  important  parameters  in  this  study  are  related  to  three 
important  controlling  parameters  :  density,  gravity  and  shear  modulus. 
Thus,  it  may  lead  us  to  three  options  to  simulate  the  gravity  effects  in 
small  scale  modeling  :  controlling  density  (p) ,  gravity  (g)  and  shear 
modulus  (G) . 

B.  RECOMMENDATIONS 

This  project  raised  the  possibility  of  further  research  topics 
relating  to  the  gravity  effects  in  small-scale  modeling. 

1.  Since  the  initial  intent  of  this  study  was  to  determine  the 
parameters  that  govern  the  gravitational  effects,  future  efforts  should 
consider  different  magnitude  of  the  each  gravity  parameter  for  more 
comprehensive  study. 
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2.  Additional  work  should  be  performed  to  investigate  the  effects 
of  gravity  parameters  in  a  model  with  a  small-scale  structure.  This 
will  help  to  determine  the  feasibility  of  controlling  gravity  parameters 
in  small-scale  modeling. 

3.  The  gravity-controlled  small-scale  model  study  without  using 
the  centrifuge  techniques  should  be  supplemented  by  a  centrifuge  model 
study  for  comparison  purpose. 

4.  A  field  study  is  necessary  to  verify  the  results  of  the 
centrifuge  model  study  and  gravity- controlled  small-scale  model  study 
once  reasonable  information  is  obtained. 
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